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predict  the  profile  using  standard  meteorological  profile  data  was  also  being 
Investigated.  The  study  features  two  modes  of  data  archiving:  (1)  continuous 
archiving  of  1  hr  average  wind  profiles  and  turbulence  levels,  and  (2)  high 
time  resolution  measurements  In  association  with  other  measurements 

(ground-based  optical  scintillometers,  aircraft  or  radiosondes).,  In  an 

intensive  'campaign'  given  the  acronym  EWAK. 

^The  atmospheric  turbulence  profiles  and  resultant  optical  propagation 

parameters  have  been  found  to  be  strongly  Influenced  by  synoptic  conditions. 
In  particular,  the  turbulence  was  substantially  affected  by  to  strength  and 
location  of  the  Jetstream.  A  very  strong  correlation  between  wind  shear 
(which  was  maximum  above  and  below  the  core  of  the  Jet)  and  pilot  reports  of 
turbulence  was  found.  Richardson  number  gave  a  much  weaker  Indication, 
possibly  because  of  the  poorer  quality  of  the  vertical  temperature  gradient 
data.^A  comparison  of  five  different  methods  (four  measurement  and  one  model) 
of  obtaining  optical  Cn2  showed  average  disagreements  as  large  as  a  factor 
of  thre^r— ^A  study  of  the  ratio  of  temperature  to  velocity  microturbulence 
showed  that  the  assumption  of  a  constant  mixing  efficiency  (used  In  the  Van 
Zandt  model)  may  not  be  valid  for  very  weak  turbulence.  ^_[he  potential  for 
using  operational  numerical  forecast  models  to  compute  turbulence  estimates 
from  predicted  wind  and  temperature  profiles  was  examined  In  a  preliminary 
look  at  the  ability  of  the  NMC  neated  Grid  Model  (NGM)  to  reproduce  the  wind 
speed  and  direction  directly  measured  by  the  radars.  The  standard  deviation 
between  the  radar  and  model  was  on  the  order  of  6  m/s  for  wind  speed  and  15 
degrees  for  wind  direction  at  the  Initial  analysis  time.  The  uncertainty  In 
wind  direction  Increased  to  about  25  to  30  degrees  for  a  48  hr  forecast  but 
the  uncertainty  In  wind  speed  did  not  change  significantly.  A  systematic 
difference  of  several  m/s  was  found  during  the  winter,  probably  due  to  over 
smoothing  of  mesoscale  features  by  the  model. 
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This  project  is  concerned  with  the  relation  of  meteorological  conditions 
to  parameters  and  processes  that  influence  the  optical  propagation  properties 
of  the  turbulent  atmosphere.  The  approach  is  centered  around  the 
establishment  of  a  climatology  of  refractive  index  structure  function 
parameter  as  measured  with  a  network  of  doppler  radars.  The  relation  of  the 
atmospheric  turbulence  profile  to  the  synoptic  context  and  physical  models  to 
predict  the  profile  using  standard  meteorological  profile  data  is  also  being 
investigated.  The  study  features  two  modes  of  data  archiving:  (1)  continuous 
archiving  of  1  hr  average  wind  profiles  and  turbulence  levels,  and  (2)  high 
time  resolution  measurements  in  association  with  other  measurements 
(ground-based  optical  scintillometers,  aircraft  or  radiosondes)  in  an 
intensive  'campaign'  given  the  acronym  EWAK. 

The  atmospheric  turbulence  profiles  and  resultant  optical  propagation 
parameters  have  been  found  to  be  strongly  influenced  by  synoptic  conditions. 

In  particular,  the  turbulence  is  substantially  affected  by  to  strength  and 
location  of  the  jetstream.  A  very  strong  correlation  between  wind  shear 
(which  is  maximum  above  and  below  the  core  of  the  jet)  and  pilot  reports  of 
turbulence  was  found.  Richardson  number  gave  a  much  weaker  indication, 
possibly  because  of  the  poorer  quality  of  the  vertical  temperature  gradient 
data.  A  comparison  of  £ive  different  methods  (four  measurement  and  one  model) 
of  obtaining  optical  C  showed  average  disagreements  as  large  as  a  factor 
of  three.  A  study  of  tne  ratio  of  temperature  to  velocity  microturbulence 
showed  that  the  assumption  of  a  constant  mixing  efficiency  (used  in  the  Van 
Zandt  model)  may  not  be  valid  for  very  weak  turbulence.  The  potential  for 
using  operational  numerical  forecast  models  to  compute  turbulence  estimates 
from  predicted  wind  and  temperature  profiles  was  examined  in  a  preliminary 
look  at  the  ability  of  the  NMC  Nested  Grid  Model  (NGM)  to  reproduce  the  wind 
speed  and  direction  directly  measured  by  the  radars.  The  standard  deviation 
between  the  radar  and  model  was  on  the  order  of  6  m/s  for  wind  speed  and  15 
degrees  for  wind  direction  at  the  initial  analysis  time.  The  uncertainty  in 
wind  direction  increased  to  about  25  to  30  degrees  for  a  48  hr  forecast  but 
the  uncertainty  in  wind  speed  did  not  change  significantly.  A  systematic 
difference  of  several  m/s  was  found  during  the  winter,  probably  due  to  over 
smoothing  of  mesoscale  features  by  the  model. 


Final  Report 


Contract  AFOSR-86-0049 

Long  Term  Studies  of  the  Refractive  Index  Structure  Parameter 
in  the  Troposphere  and  Stratosphere 

C.W.  Fairall,  D.W.  Thomson,  and  W.J.  Syrett 
Department  of  Meteorology 
505  Walker  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Aooesslon  For 

_ l _ 

NTIS  GRA&I  fij 

DTIC  TAB  U 

Unannounced 

□ 

Justification — 

By - - 

Dlstrlbut  Ion/ _ _ 

Availability  Codes 
(Avai'l  and/or 


Dist 


Special 


This  project  is  concerned  with  the  relation  of  meteorological  conditions 
to  parameters  and  processes  that  influence  the  optical  propagation  properties 
of  the  turbulent  atmosphere.  The  approach  is  centered  around  the 
establishment  of  a  climatology  of  refractive  index  structure  function 
parameter  as  measured  with  a  network  of  Doppler  radars .  The  relation  of  the 
atmospheric  turbulence  profile  to  the  synoptic  context  and  the  use  of  physical 
models  to  predict  the  profile  using  standard  meteorological  profile  data  is 
also  being  investigated.  The  study  features  two  modes  of  data  archiving:  (1) 
continuous  archiving  of  1  hr  average  wind  profiles  and  turbulence  levels,  and 
(2)  high  time  resolution  measurements  in  association  with  other  measurements 
(ground-based  optical  scintillometers,  aircraft  or  radiosondes)  in  an 
intensive  'campaign'  given  the  acronym  EWAK. 

The  atmospheric  turbulence  profiles  and  resultant  optical  propagation 
parameters  have  been  found  to  be  strongly  influenced  by  synoptic  conditions. 

In  particular,  the  turbulence  is  substantially  affected  by  to  strength  and 
location  of  the  jet  stream.  A  very  strong  correlation  between  wind  shear 
(which  is  maximum  above  and  below  the  core  of  the  jet)  and  pilot  reports  of 
turbulence  was  found.  Richardson  number  gave  a  much  weaker  indication, 
possibly  because  of  the  poorer  quality  of  the  vertical  temperature  gradient 
data.  A  comparison  of  five  different  methods  (four  measurement  and  one  model) 
of  obtaining  optical  Cn  showed  average  disagreements  as  large  as  a  factor 
of  three.  A  study  of  the  ratio  of  temperature  to  velocity  microturbulence 
showed  that  the  assumption  of  a  constant  mixing  efficiency  (used  in  the  Van 
Zandt  model)  may  not  be  valid  for  very  weak  turbulence.  The  potential  for 
using  operational  numerical  forecast  models  to  compute  turbulence  estimates 
from  predicted  wind  and  temperature  profiles  was  examined  in  a  preliminary 


look  at  the  ability  of  the  NMC  Nested  Grid  Model  (NGM)  to  reproduce  the  wind 
speed  and  direction  directly  measured  by  the  radars .  The  standard  deviation 
between  the  radar  and  model  was  on  the  order  of  6  m/s  for  wind  speed  and  15 
degrees  for  wind  direction  at  the  initial  analysis  time.  The  uncertainty  In 
wind  direction  increased  to  about  25  to  30  degrees  for  a  48  hr  forecast  but 
the  uncertainty  in  wind  speed  did  not  change  significantly.  A  systematic 
difference  of  several  m/s  was  found  during  the  winter,  probably  due  to  over 
smoothing  of  mesoscale  features  by  the  model. 
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[ntroduction 


The  importance  of  turbulence  in  the  atmospheric  boundary  layer,  which  is 


the  basic  living  environment  of  humankind,  has  long  been  recognized.  Only  in 


the  last  few  decades  has  the  importance  of  free  atmospheric  (i.e.,  above  the 


boundary  layer)  turbulence  been  recognized.  Even  in  the  absence  of  clouds, 


stably  stratified  fluids  (e.g.,  the  free  troposphere  and  stratosphere)  are 


observed  to  experience  intermittent  transitions  to  turbulent  regimes. 


Although  it  is  referred  to  as  Clear  Air  Turbulence  (CAT)  it  does  occur  in 


regions  of  stable  clouds  such  alto  or  cirmstratus.  The  beautiful  breaking 


wave  clouds  provide  visual  evidence  of  CAT.  CAT  is  important  in  aircraft 


performance  and  safety,  optical  propagation,  EM  propagation  and 


communications,  atmospheric  dispersion  of  pollutants,  and  as  a  source  of 


viscous  drag  in  the  'free'  atmosphere. 


The  atmospheric  refractive  index  structure  function  parameter,  C 


is  important  in  a  broad  range  of  optical  propagation  applications  (see  Section 


II) .  Examples  of  the  consequence  of  atmospheric  microturbulence  on  optical 


systems  include:  reduced  intensity  (Yura,  1971),  beam  wander  (Fried,  1966), 


scintillation  and  coherence  (Fried  and  Schmeltzer , 1967) ,  and  anisoplanatism 


(Fried,  1981) .  The  consequences  of  spatially  and  temporally  varying  C 


can  be  computed  for  a  specific  optical  system,  at  least  in  principle,  by 


specifying  the  values  of  Cn  along  the  optical  path.  While  an  extensive 


base  of  data  and  theory  is  available  for  planetary  boundary  layer  (PBL) 


properties  of  Cn  (see  Fairall  et  al.,  1982;  Fairall,  1987),  only  a  few 


case  studies  are  available  for  the  free  troposphere  and  the  lower  stratosphere 


(Walters  and  Kunkel,  1981).  The  ALLCAT  (i.e.,  HICAT,  MEDCAT,  etc.;  program  of 


the  late  1960 's  focused  primarily  on  large  scale  turbulence  severe  enough  to 


damage  aircraft  or  discomfit  passengers.  Now  the  structure  and  dynamical 
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properties  of  both  CAT  and  in  the  free  troposphere  and  lower 

2 

stratosphere  are  receiving  renewed  attention.  For  this  is  because  of 

the  recognition  of  the  importance  of  the  isoplanatic  angle  in  the  performance 

of  a  variety  of  ground- to- space  systems  that  utilize  adaptive  optics. 

Concurrently  development  of  VHF  clear  air  Doppler  radars ,  the  performance  of 

2 

which  are  directly  proportional  to  Cn  ,  has  rekindled  the  interest  of  the 

meteorological  community  in  the  climatology  and  dynamics  of  CAT.  It  turns  out 

that  this  is  still  poorly  understood.  For  example,  a  recent  long  term  study 
2 

of  Cn  using  VHF  radar  (Nastrom  et  al.,  1981)  revealed  substantial  diurnal 
and  seasonal  variations  even  in  the  stratosphere.  This  study  suggested  that 
strong  tropospheric  convection  and  the  jet  stream  were  relevant  factors  in  the 
intensity  of  CAT  but,  at  present,  this  is  only  conjecture  based  on  physical 
plausibility. 

The  investigation  of  atmospheric  turbulence  involves  the  collection  of 

data,  the  development  of  theory  and  the  implementation  of  models.  The  data 

tells  us  what  is  'up  there'  at  the  time  of  the  measurement,  but  the  theory  and 

models  are  true  indices  of  our  understanding  of  the  phenomena.  With  models  we 

attempt  to  predict  an  important  but  difficult  to  measure  variable  (e.g., 

2 

Cn  )  from  riables  that  we  expect  to  have  at  our  disposal  (e.g.,  regional 
scale  radiosonde  information) .  Dynamical  models  of  turbulence  based  on  the 
Navier-Stokes  and  turbulent  kinetic  energy  (TKE)  equations  have  a  rich 
tradition  in  PBL  research.  Both  second  order  closure  and  large  eddy 
simulations  have  provided  great  insights  into  the  structure  of  PBL  turbulence. 
In  contrast,  the  only  model  of  CAT  in  use  today  is  purely  statistical  in 
nature  (Van  Zandt  et  al. ,  1981;  Moss,  1986). 

Remote  sensors  are  ideal  for  the  study  of  climatological  properties  of 
turbulence  in  the  free  troposphere  and  stratosphere.  In  this  regard,  the  VHF 
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and  UHF  Doppler  radars  are  particularly  well  suited  because  of  their  all 
weather  and  day /night  capabilities.  The  Penn  State  Department  of  Meteorology 
is  now  completing  a  mesoscale  triangular  network  of  three  VHF  radars  under 
funding  provided  by  the  DoD  (AFOSR)  University  Research  Instrumentation 
Program.  DoD  (AFOSR  and  ONR)  funding  has  also  been  obtained  for  the 
construction  of  a  multichannel  mm-wave  radiometer  system  for  continuous 
groundbased  measurements  of  temperature  profiles  and  integrated  water 
vapor/liquid.  Subsystems  of  this  radiometer  are  now  being  tested.  It  is 
expected  to  be  online  at  Penn  State  in  mid  1988.  Finally,  a  fourth  Doppler 
radar  has  been  constructed.  It  is  a  UHF  system  similar  in  concept  to  the  VHF 
systems  but  operating  at  0.75  m  rather  than  6  m  wavelength.  The  UHF  system 
has  a  smaller,  more  portable  antenna,  better  vertical  resolution  (100  m  vs. 

300  m  for  the  VHF  systems)  and  far  superior  low  altitude  capabilities  (minimum 
range  of  200  m) .  The  radar  systems  are  capable  of  producing  a  wind  and 
turbulence  profile  with  roughly  30  to  90  second  time  resolution.  Thus,  they 
are  also  ideal  for  highly  detailed  studies  of  turbulence  and  its  relationship 
to  mesoscale  phenomena.  The  combination  of  the  mesoscale  triangle  of  VHF 
profilers,  the  UHF  profiler  and  the  radiometers  constitutes  an  atmospheric 
observing  system  that  represents  a  quantum  jump  in  our  ability  to  study  the 
atmosphere . 

This  report  describes  a  two  year  project  to  use  these  systems  to  study 

2 

the  seasonal  and  diurnal  climatology  of  ,  the  atmospheric  dynamical 

2 

processes  responsible  for  the  variability  of  Cn  ,  and  to  investigate 

2 

measurement  methods  and  models  of  Cn  .  The  results  have  been  published  in 
the  open  literature  (see  section  VIII)  and  in  a  number  of  graduate  assistant 
M.S.  theses  [Beecher  (1987),  Carlson  (1987),  Knowlton  (1987),  Neiman  (1987), 
and  Syrett  (1987)].  The  work  combined  the  full  power  of  the  department's 
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observing  systems,  access  to  national  data  networks  and  weather/satellite 
information,  and  cooperation  with  a  number  of  government  research 
laboratories .  During  the  study  an  intensive  field  program  was  held  at  Penn 
State  involving  in  situ  measurements  with  aircraft  and  thermosonde  balloons 
and  a  variety  of  ground  based  electro-optical  systems  (e.g.,  scintillometers). 
This  study  was  given  the  acronym  EWAK  and  involved  scientists  from  a  number  of 
DoD  laboratories  (see  section  IV) .  The  purpose  of  EWAK  was  to  compare 
measurement  methods  and  to  test/develop  models  of  clear  air  turbulence.  The 
emphasis  on  relating  turbulence  to  mesoscale  and  synoptic  scale  structures  was 
intended  to  promote  the  extension  of  the  results  to  different  climatological 
regions .  A  thorough  understanding  of  these  relationships  should  lead  to 
assessment  and  forecasting  of  atmospheric  turbulence  with  a  combination  of 
satellite  informational  and  numerical  global  weather  models. 


For  many  purposes  the  effects  of  atmospheric  gas  on  propagating 
electromagnetic  and  acoustic  radiation  may  be  conveniently  subdivided  into  the 
following  subcategories : 

a)  mean  density  gradients  result  in  beam  refraction  and,  consequently, 
tracking  or  pointing  errors, 

b)  density  and  velocity  fluctuation  caused  by  turbulence  produce 
refractive  index  variations  which,  in  turn,  degrade  system  performance. 

Regarding  the  bulk  radar  refractive  index,  n,  of  atmospheric  gas  it  is 
convenient  to  specify  it  in  terms  of  the  refractivity ,  N  -(n-l)*10B,  in 
which  case 

N-  77.6  (P  +  4810  e/T)/T  (1) 

where  N,  X,  P  and  e  are  in  units  of  ppm,  K,  mb  (total  pressure),  and  mb  (vapor 

pressure),  respectively.  The  second  term  including  e  specifies  the 

contribution  by  polar  molecules  (principally  water).  It  is  often  negligible 

at  optical  frequencies  where  the  humidity  coefficient  is  much  smaller. 

For  specification  of  the  effects  of  turbulence  the  relevant  atmospheric 

2 

properties  are  the  refractive  index  structure  function  parameter,  Cn  , 
(Tatarski,  1961)  and  the  inner  scale  or  Kolmogorov  scale  of  turbulence, 
(Livingston,  1972;  Hill  and  Clifford,  1978).  C  2  can  be  related  to 
microraeteorological  variables  by  (Wesely,  1976) 

Cn  “  A(GT  +  2  B  CTq  +  fi2cq)  (2) 

where  A  and  B  are  functions  of  temperature  (T) ,  pressure  and  specific  humidity 


2  2 

(q) ;  C^,  ,  and  are  the  temperature  and  humidity  structure 

function  parameters.  The  values  of  A  and  B  are  also  radiation- type  and 
wavelength  dependent.  is  related  to  the  viscosity  and  density  of  air 

(both  functions  of  altitude)  and  the  rate  of  dissipation  of  turbulent  kinetic 
energy,  «  (Hinze,  1975). 


Ill .  Background  on  Clear  Air  Turbulence 


A.  Micro turbulence 

In  general,  atmospheric  turbulence  is  anisotropic.  However,  it  is  known 

that  at  small  size  scales  (large  values  of  wavenumber,  k) ,  the  eddies  become 

increasingly  isotropic.  In  the  isotropic  limit  it  can  be  shown  that  (Hinze, 

2 

1975)  the  structure  function  parameter,  Cx  ,  for  the  unspecified 
variable, X,  (X  could  be  u,  T,  q,  or  n) 

cl  -  <(X(r) -X(r+d) )2>/d2/3  (3) 

is  independent  of  the  spacing,  d.  In  this  equation  X(r)  denotes  the  value,  of 
X  at  position  r  while  X(r+d)  denotes  the  value  at  a  position  a  distance  d  from 
r;  the  brackets  denote  an  average.  Using  dimensional  arguments,  Kolmogorov 
showed  that  in  the  isotropic  limit  the  one -dimensional  variance  spectrum  obeys 
a  k  wavenumber  dependence 

*x<k>  “  0xXx«  l//3  k‘5/3  (4) 

where  xx  is  the  rate  of  dissipation  of  one  half  of  the  variance  of  X  and 
is  an  empirical  constant  (on  the  order  of  0.5  for  velocity  and  0.8  for 
scalars) .  It  can  similarly  be  shown  that 

<px(k)  -  0.25  Cx  k'5/3  (5) 

where  the  factor  0.25  represents  several  mathematical  constants.  Note  that 
(4)  and  (5)  imply  the  Corrsin  relation 
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C2  -  4  Mx^1/3 


(6) 


Velocity  is  a  special  case  (since  Xu”0  so  that 


C2  -  2.0  e2/3 
u 


(7) 


B.  TKE  and  Variance  Budget  Equations 


The  simplified,  horizontally  homogeneous  budget  equations  for  TKE  (symbol 


E)  and  scalar  variance  (symbol  <x>) ,  are 


DE/Dt  -  -<uw>3U/3z  +  (g/0)30/3z  -e  +  transport 


(8a) 


D(<x  >/2)/Dt  -  -<xw>3X/3z  -*x  +  transport 


(8b) 


where  the  small  symbols  denote  turbulent  fluctuations  and  the  capital  symbols 
denote  average  values,  U  is  the  streamwise  velocity,  W  the  vertical  velocity, 
g  the  acceleration  of  gravity,  6  the  potential  temperature  and  z  the  altitude. 
Assuming  a  state  of  dynamic  equilibrium,  neglecting  transport,  and  using  an 
eddy  diffusion  coefficient,  K,  to  express  the  covariance  term 


<uw>  -  -K  3U/3z 
m 


(9a) 
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€  -  Km(3U/3zr  [1  -  Ri/Pr] 


Xx  “  Kx(3X/3z)‘ 


(10a) 


(10b) 


where  Pr-K^/K^  is  the  turbulent  Prandtl  number  and  Ri  is  the  gradient 


Richardson  number 


Ri  -  tr/s 


The  factor  S  is  the  square  of  the  vertical  gradient  of  the  vector  mean  wind 
and  N  is  the  Brunt -Vaisala  frequency,  which  in  the  atmosphere  is  approximately 
given  by 


*T  -  ( g/9 )  30/3z 


Using  (10)  we  can  obtain  the  ratio  of  the  scalar  dissipation  to  the  TKE 


dissipation 


*x/€  “  (W  Cx  /Cu  "  [  (3X/3z)  2/N2  ]  Ri/(Pr-Ri) 


In  the  oceanographic  literature  (Gregg,  1987),  this  is  expressed  in  terms  of 


the  mixing  coefficient,  7, 


7X  ”  XxN2/[(aX/3z)2e] 


In  actively  turbulent  layers  with  Pr=*l  and  Ri=»0.25,  (13)  implies 
7-Ri/(Pr-Ri)=»l/3 .  Gossard  and  Frisch  (1987)  have  used  the  shear  budget 


equations  to  show  that  ^x7//3u=*3/2  while  Gage  and  Nastrom  (1986)  have  shown 
that  equipartition  of  two-dimensional  TKE  and  potential  energy  due  to  gravity 
wave  displacements  in  a  stratified  fluid  implies  for  large 

scale  (anisotropic)  turbulence.  As  we  shall  discuss  later,  relations  derived 
from  (13)  have  been  used  to  compute  e  and  Km  from  clear  air  radar  data. 

C.  Clear  Air  Turbulence  Length  Scales 

The  fif'd  of  turbulence  contains  a  bewildering  variety  of  length  scales 
which  we  will  not  attempt  to  discuss  here.  For  our  purposes,  we  focus  on  the 
smallest  scales  present,  the  Kolmogorov  microscale  L^,  and  the  scale  of  the 
dominant  vertical  motions  (eddy  overturning  scale)  Lq,  which  is  related  to 
the  integral  scale  calculated  from  the  vertical  velocity  autocorrelation 
function.  At  the  Kolmogorov  scale,  viscosity  is  rapidly  destroying  the 
turbulent  fluctuations;  the  spectrum  begins  to  deviate  from  the  -5/3  behavior 
at  sizes  an  order  of  magnitude  larger  than  this  (Hill  and  Clifford,  1978) .  In 
stratified  turbulence,  the  overturning  scale  is  proportional  to  the  buoyancy 
length  scale,  (Gregg,  1987),  also  referred  to  as  the  Ozmidov  scale.  In 
terms  of  the  turbulent  parameters,  these  scales  are 

1^  -  (u/«3)1/4  (15a) 

1^  -  (£/N3)1/2  (15b) 

The  inertial  subrange  occurs  for  those  size  scales  smaller  than  and 
greater  than  L^.  If  L^L^,  then  the  energy  containing  vertical  motions 
are  rapidly  destroyed  by  viscosity;  if  L^^IOL^,  then  no  inertial  subrange 
is  anticipated. 

The  ratio  of  the  length  scales  forms  a  natural  activity  parameter  to 


L2- 


classify  the  turbulent  state  (Gibson,  1987) 


A  -  ( W 


V3  .  £/(uN2) 


Gregg  (1987)  suggests  the  following  empirical  classifications 


Value  of  A 


Turbulent  State 


Decaying  turbulence ,  <wu>— <w0>=*O 


A>200 


Isotropic 


A>10000 


Fully  developed 


The  physical  interpretation  of  the  mixing  coefficient  and  the  activity 
parameter  can  be  illustrated  by  noting  that  the  turbulent  diffusivity  can  be 
expressed  as  K^-A^^u.  Thus,  when  the  product  of  the  activity  parameter 
and  the  mixing  coefficient  exceeds  one,  then  the  turbulent  mixing  processes 
are  more  efficient  than  molecular  diffusivity. 

The  concept  of  length  scales  is  also  used  to  eliminate  the  eddy  diffusion 
coefficients  from  (10)  by  invoking  mixing  length  theory  (Hinze,  1975) 


K  -  c,E1/2  L 
1  o 


where  c^  is  an  empirical  constant  and  the  square  root  of  E  represents  a 
velocity  scale.  The  dissipation  is  also  related  to  these  quantities 


‘  -  c~E  '  /h 


These  relations  can  be  used  in  (10)  to  yield 


Vi? 
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-  2  *2/3-  2  (c-j/c^3)  (3U/3z)2  L^/3 

-  2  (£x/£u)  (c1/c2/3^/ (Pr-Ri)  (ax/az)2 


Notice  that  (19a)  can  be  manipulated  to  give 


e/N 
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(c3/2/cV2)/Ri3/2  L2 


(19a) 

(19b) 

(20) 
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The  most  common  convention  is  to  fix  the  ratio  c-^/Cj  '  -1  and  to  assume 
that  in  actively  turbulent  layers  Ri“0.25.  This  implies  that  Lq-0.35 
(Hocking,  1982).  Because  the  constants  are  chosen  arbitrarily,  there  is  no 
physical  significance  to  this  particular  ratio  of  Lq  to  1^.  Since  is 
considered  to  be  the  outer  limit  of  the  validity  of  the  inertial  subrange,  it 
makes  sense  that  the  true  integral  scale  and  the  energy  containing  scale  are 
larger  than  L.  .  It  is  also  clear  that  application  of  these  expressions  is 
likely  to  be  confused  by  our  inability  to  be  sure  of  the  values  of  Pr  and  Ri 
and  by  the  fact  that  the  assumptions  of  stationarity  and  negligible  transport 

will  not  be  valid  in  all  conditions. 

2 

D.  Statistical  Models  of  C 

n 

2 

Van  Zandt  et  al.  (1978,1981)  have  developed  a  model  of  Cn  based  on 
a  statistical  integration  of  simplified  forms  of  (19).  In  this  approach,  a 
smooth  mean  profile  with  a  vertical  resolution  roughly  equivalent  to  a 
rawinsonde  is  the  input.  Velocity  and  temperature  (and,  therefore,  shear  and 
temperature  gradient)  fluctuations  are  allowed  (relative  to  the  mean  profile) 
with  probability  distributions  obtained  from  a  mix  of  empirical  analysis  and 
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theory  about  gravity  wave  effects.  Active  turbulent  regions  are  assumed  to 


Si 
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exist  when  the  fluctuations  produce  Ri<0.25.  By  integrating  the  joint 

probability  distribution  over  shear,  temperature  gradient,  and  size  scale 

2 

space,  an  average  value,  <Cn  >,  is  computed. 

2 

The  value  of  <Cn  >  obtained  in  this  manner  is  the  expected  value  at 

some  altitude;  as  such,  it  does  not  actually  contain  any  information  about  the 

vertical  distribution.  We  can  see  that,  at  any  specific  time,  quite  different 
2 

results  for  Cn  can  be  obtained  from  measurements  depending  on  the 
vertical  resolution.  Suppose  that  the  vertical  resolution  (e.g.,  a  radar 
range  gate  of  several  hundred  meters)  is  much  greater  than  the  turbulent  patch 
thickness,  H.  In  this  case,  we  would  expect  that  the  radar  produces 
sufficient  averaging  to  be  consistent  with  model  assumptions.  A  high 
resolution  radar,  or  a  high  vertical  resolution  in  situ  measurement  (e.g.,  an 
aircraft  flying  at  constant  altitude)  is  likely  to  produce  a  measurement  that 
is  either  in  an  inactive  region  or  in  an  active  region.  Thus  a  high  resolution 
measurement  at  a  particular  altitude  may  require  averaging  over  a  long  period 
of  time  to  be  consistent  with  the  model.  The  time  required  would  be  many 
times  longer  than  the  typical  lifetime  of  an  actively  turbulent  layer,  which 
appears  to  be  on  the  order  of  a  few  hours  (Syrett,  1987;  also,  see  the 
discussion  in  Section  III-F). 

In  an  earlier  paper  Van  Zandt  et  al.  (1978)  examined  this  issue  with 
their  statistical  model  by  computing  the  probability  of  turbulence  occurring 
at  a  fixed  altitude,  F.  This  has  been  interpreted  by  some  (e.g.,  Weinstock, 
1981)  as  the  expected  fraction  of  a  range  gate  of  thickness,  2A,  that  is 
occupied  by  turbulence  of  average  scale,  L^“Lq,  i.e.,  F=‘L^/(2A) . 

Obviously,  this  assumes  that  the  patch  thickness  is  much  less  than  the  range 
gate  thickness  and  that  only  one  layer  is  likely  to  occur  per  range  gate. 

Thus,  average  measurements  produced  in  this  fashion  are  interpreted  as 


<€>  -  F  € 


(21a) 


<C^>  -  F  C*  (21b) 

where  the  values  on  the  right  hand  side  represent  data  within  the  active 

layers  where  Ri“0.25.  While  the  value  of  F  depends  on  the  mean  conditions, 

1/3  2 

Gage  et  al.  (1980)  use  the  model  to  show  that  F  '  N  is  approximately 
constant  and  has  the  value  of  4*10 in  the  troposphere  and  8*10"^  in  the 
stratosphere. 

The  Van  Zandt  model  has  never  been  tested  in  detail.  It  has  been 

2 

evaluated  by  comparing  radar  measured  with  predictions  from 

rawinsondes.  In  a  few  cases,  optical  and  in  situ  data  have  also  been  used. 

2 

On  average,  the  model  does  quite  well  for  Cn  ;  this  is  not  shocking 

because  the  model  originally  contained  one  unspecified  constant  which  was 

selected  to  fit  radar  data.  The  internal  details  of  the  model,  such  as  the 

probability  distributions  for  shear,  N,  and  F,  have  not  been  evaluated.  Also, 

the  model  appears  to  overpredict  e  by  two  orders  of  magnitude  (Fairall  and 

Markson,  1984) .  Even  the  model  predictions  of  have  never  been 

examined  for  a  large  data  set.  Furthermore,  the  model  is  incomplete  in  that 

it  provides  information  about  the  probability  of  observing  turbulence  (F)  and, 

perhaps,  the  typical  size  scale  of  the  turbulence  (<Lq>) ,  but  it  tells  us 

nothing  about  the  vertical  distribution  of  turbulent  regions  nor  does  it 

guarantee  that  the  patch  size  (H)  is  the  same  as  Lq.  Observations  in  the 

ocean  (Gregg,  private  communication)  and  in  the  atmosphere  (Barat  and  Bertin, 

1984)  show  that  H»L  . 

o 


E.  Inactive  Regions 


The  Van  Zandt  model  partitions  the  atmosphere  vertically  into  regions 


with  active  turbulence  (Ri<0.25)  and  regions  that  are  considered  to  be 

nonturbulent  (or  the  turbulence  makes  a  negligible  contribution  to  the  average 
2 

of  C  ) .  By  inactive,  we  mean  that  there  is  no  production  of  TKE  or 
variance  which  implies  that  the  covariance  terms  in  (9)  must  be  zero.  This 
condition  is  met  when  the  activity  parameter  is  less  than  15.  Thus,  we  can 
define  a  threshhold  value  (Gibson,  1987)  for  the  dissipation  rate,  et,  such 
that  values  of  e  less  than  e t  imply  decaying  turbulence, 


et/(utO  -  15 


-  8  2-3 

A  typical  value  in  the  free  troposphere  is  et=«5*10’  ms  .  Even  if 
A>15  the  turbulence  may  be  decaying  if  the  destruction  terms  exceed  the 
production  terms.  Sometimes  the  term  'fossil'  turbulence  is  used  to  describe 
this  state  although  that  term  is  also  applied  to  residual  temperature 
structure  that  remains  after  the  velocity  turbulence  has  been  consumed. 

Decaying  turbulence  has  been  extensively  studied  in  the  laboratory  using 
flow  through  grids  (e.g.,  Itsweire  et  al. ,  1986).  Atmospheric  studies  are 
quite  rare,  although  Nieuwstadt  and  Brost  (1986)  performed  a  large  eddy 
simulation  model  study  of  the  decay  of  convective  turbulence.  We  can  use  the 
TKE  budget  equation  to  analyze  decaying  turbulence,  where  we  drop  the 
production  terms 


dE/dt  -  -€ 


However,  if  we  eliminate  E  with  (18),  then  we  obtain  a  differential  equation 
for  e 
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de/e4/3  -  -(3/2)  (c2/Lq)2/3  dt 


If  we  assume  that  the  size  scale  remains  constant  throughout  the  decay  process 
then  (24)  has  the  solution 


€  -  €q  (1  +  t/tQ) 


where  e  is  the  value  of  e  at  the  beginning  of  the  decay  process  and  t  is 


„  0  rT2.,  2  ,,1/3  -2/3  „-l 

“  2  tV(c2  €o)]  c?  N 


Measurements  in  the  ocean  (Dillon,  1982)  suggest  tQN=«0.3,  which  implies  that 
an  inactive  region  has  a  turbulent  lifetime  that  is  only  on  the  order  of  a 
Brunt -Vaisala  period.  However,  this  is  the  value  of  Brunt -Vaisala  period 
within  the  turbulent  layer  itself,  which  may  be  much  longer  (because  of 
mixing)  than  the  value  associated  with  the  mean  background  temperature  profile 
(about  two  minutes  in  the  atmosphere) . 

F.  Active  Regions 

The  dynamics  and  time  scales  of  regions  of  active  turbulence  are  fairly 
complicated  and  very  little  complete  data  is  available  to  aid  in  the  analysis. 
It  is  generally  conceded  that  the  turbulence  is  initiated  by  Kelvin -Helmholtz 
instability  when  the  Richardson  number  decreases  to  a  value  near  0.25. 

Several  physical  processes  are  able  to  cause  this  decrease.  Nonturbulent 
portions  of  the  atmosphere  are  subjected  to  fluctuations  in  shear  and/or 
potential  temperature  gradient  caused  by  gravity  wave  disturbances.  These 
fluctuations  may  cause  turbulence  as  described  in  the  section  on  statistical 
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models  of  CAT.  Synoptic  and  mesoscale  dynamical  processes  also  cause 
evolutions  in  the  shear  and  lapse  rates  which  can  cause  turbulence.  A  classic 
example  of  this  is  the  persistent  regions  of  CAT  found  above  and  below  well 
developed  jet  streams.  If  we  take  the  vertical  derivatives  of  the  standard 
budget  equations  for  the  mean  wind  and  potential  temperature,  then  we  can 
create  budget  equations  for.  the  evolution  of  mean  shear  and  mean  lapse  rate. 
These  can,  in  principle,  be  combined  to  yield  a  budget  equation  for  the 


Richardson  number-  a  ' Richardsonnuraberogenesis '  process,  to  paraphrase 
meteorological  jargon. 

Once  a  layer  of  thickness,  H,  becomes  turbulent,  it  is  of  interest  to 
ponder  the  temporal  duration  of  the  turbulent  event.  If  the  turbulence  is 
caused  by  a  gravity  wave  fluctuation,  then  we  do  not  expect  the  event  to  last 
more  than  a  fraction  of  an  inertial  period  (say,  a  few  hours).  If  the 
breakdown  is  due  to  synoptic  processes,  then  in  principle  the  turbulence  can 
endure  as  long  as  the  'Richardsonnumberogenesis'  can  maintain  the  instability 
against  turbulent  mixing.  Once  a  layer  is  turbulent,  the  mixing  process  tends 
to  cause  Ri  to  increase.  This  is  because  3U/3z  an  dd/dz  are  both  reduced  by 
mixing,  but  3U/3z  occurs  as  the  inverse  square  in  Ri.  Thus,  we  anticipate 
that  the  mixing  caused  by  the  turbulence  forces  Ri  to  increase  until  Ri=l.  At 
this  point,  shear  production  of  TKE  is  roughly  canceled  by  buoyant  destruction 
and  the  turbulence  begins  to  decay  due  to  dissipation.  This  process  occurs 
long  before  we  reach  et>  In  other  words,  shear  production  is  still  active, 
the  fluxes  are  still  nonzero,  but  buoyant  destruction  and  dissipation  exceed 
production.  The  layer  is  decaying  but  it  is  still  active. 

This  suggests  that  we  must  consider  two  more  time  scales:  the  time  scale 
for  the  mixing  process  to  occur,  r  ,  and  the  time  scale  for  the  active 
region  to  decay  to  We  assume  that  the  mixing  process  must  transform 


kinetic  energy  into  potential  energy  by  destroying  the  ambient  potential 

temperature  gradient  for  a  layer  of  thickness.  Gregg  (1987)  has  shown  that 

2  2 

this  requires  an  amount  to  energy  equal  to  (1/12)N  H  ,  where  N  is  computed 
from  the  background  lapse  rate  when  turbulent  breakdown  occurs .  Suppose  we 
let  P  be  the  shear  production  of  TKE  and  B  the  buoyant  destruction  integrated 
of  the  entire  thickness  of  the  layer  (thus,  transport  terms  become  zero) ;  then 
in  a  state  of  slowly  evolving  equilibrium,  the  TKE  budget  equation  implies  the 
balance 


P  -  B  +  e 


If  we  substitute  the  flux  Richardson  number,  Rf-P/B,  we  get  a  relation  for  the 
mixing  time  constant 


erm  -  (l/Rf-l)Brm  -  (1/Rf- 1)N2H2/12  -  [(Pr-Ri)/Ri]  N2H2/12  (28) 


Therefore,  we  can  write  r  as 

m 


rm-  N2H2/ (12  «  7) 


Aircraft  measurements  of  the  covariance  and  gradient  terms  (Kennedy  and 
Shapiro;  1975,  1980)  in  the  vicinity  of  the  jet  stream  give  an  average  value 
of  7«0.40.  Representative  values  in  (29)  give  on  the  order  of  1  hour  for 
a  100  m  thick  layer  with  e-l*10  ^  m2s  ^ .  Notice  that  (29)  implies  that 
thicker  layers  and  weaker  layers  («  smaller)  will  persist  much  longer. 

Once  the  layer  has  been  mixed  as  discussed  above,  then  the  turbulence 
begins  to  decay  with  the  time  constant  t  described  in  III-D.  According  to 
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The  use  of  radar  data  in  optical  propagation  and  CAT  studies  requires 
consideration  of  several  factors:  radar  resolution,  range  and  measured 
variables,  the  relationship  ot  radar  and  optical  parameters,  and  the  types  of 
studies  that  are  appropriate  for  these  systems.  Details  are  provided  in  the 
literature  (e.g.,  Hocking,  1982)  and  in  our  previous  proposal,  but  a  brief 
summary  of  the  main  points  will  be  presented  here. 

The  50  MHz  radars  have  vertical  range  gate  resolution  of  300  m  up  to  8  km 
and  900  m  resolution  up  to  18  km.  The  high  resolution  mode  requires  90  s  to 
obtain  one  profile,  the  low  resolution  3  minutes.  In  normal  operations,  a 
pair  of  resolution  profiles  is  obtained  in  5  minutes.  The  UHF  radar  has  100  m 
resolution  to  2.6  km,  300  m  resolution  to  8  km,  and  900  m  resolution  to  12  km. 
Both  systems  have  two  horizontal  and  one  vertical  beam. 

The  basic  raw  data  produced  by  the  radar  is  the  mean  Doppler  shift  (which 
is  used  to  compute  the  mean  wind  vector) ,  the  width  of  the  Doppler  spectrum 
(which  is  related  to  the  turbulence  associated  with  the  wind  variance  or  the 
mean  shear  and  can  be  used  to  estimate  «),  and  the  backscatter  intensity  of 

the  signal  (which  is  related  to  the  radar  cross  section  and  is  used  to 

2  2 
calculate  radar  Cn  ) .  Details  concerning  the  computation  of  Cn  from 

the  radar  signal  are  discussed  in  Appendix  B.  Thus  the  radar  provides 

quantitative  measurements  of  the  profiles  of  U,  3U/3z,  a  ,  ow,  and 
2 

Cn  .  Models  can  be  applied  to  this  data  to  estimate  e  (e.g.,  equation 

13),  K  (e.g.,  equation  10),  and  even  d$/dz. 

The  radars  normally  operate  continuously  and  archive  data  in  one  hour 

average  blocks.  For  intensive  experiments,  we  archive  the  high  resolution 

(temporally  and  spatially)  data,  including  the  raw  Doppler  spectra  and  all 

2 

moment  tables.  An  example  of  a  time  series  of  Cn  from  two  range  gates  of 


a  VHF  system  is  shown  in  Fig.  1.  A  few  examples  of  studies  ideally  suited  to 

2 

these  systems  are:  vertical  distribution  of  Cn  'hot  spots' ,  comparison  of 

2 

daytime  versus  nighttime  Cn  at  different  times  of  the  year,  coincidence 
of  high  turbulence  levels  and  wind  shear/jet  stream  regions,  and  correlation 
with  synoptic  regimes. 
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A .  Background 

This  project  is  concerned  with  the  relationship  of  meteorological 


conditions  to  parameters  and  processes  that  influence  the  optical  propagation 
properties  of  the  turbulent  atmosphere .  The  approach  is  centered  around  the 
establishment  of  a  climatology  of  refractive  index  structure  function 
parameter  as  measured  with  a  network  of  Doppler  radars .  The  relation  of  the 
atmospheric  turbulence  profile  to  the  synoptic  context  and  a  model  to  predict 
the  profile  using  standard  meteorological  profile  data  is  also  being 
investigated.  The  study  features  two  modes  of  data  archiving:  (1)  continuous 
archiving  of  1  hr  average  wind  profiles  and  turbulence  levels,  and  (2)  high 
time  resolution  measurements  in  association  with  other  measurements 
(ground-based  optical  scintillometers,  aircraft  or  radiosondes)  in  an 
intensive  'campaign'  given  the  acronym  EWAK. 

This  project  funded  one  graduate  assistant  (Syrett)  who  studied  the 
correlation  of  turbulence  with  synoptic  context  (in  relation  to  jet  streams). 
Two  Air  Force  graduate  students  have  also  been  working  on  the  project.  Capt. 
Michael  Moss  is  completing  a  Ph.D.  dissertation  involving  radar  data  and  the 
Van  Zandt  model  and  2nd  Lieut.  Elizabeth  Beecher  has  completed  a  comparison 
(M.S.  thesis)  of  meteorological  and  optical  microstructure  from  radar  and  in 
situ  (particularly  aircraft)  measurements. 

B.  The  EWAK  Experiment 

A  major  optical/meteorological  experiment  (acronamed  EWAK)  was  held  at 
Penn  State  from  15  April  to  15  May,  1986.  Rome  Air  Development  Center  (RADC) 
was  planning  a  combined  aircraft/optical  experiment  for  the  spring  and  we 
convinced  them  to  hold  the  field  program  at  PSU  in  order  to  take  advantage  of 
the  radar  data.  Given  the  scale  of  this  experiment,  optical  propagation 
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scientists  from  several  other  laboratories  were  also  invited  to  participate. 
The  enclosed  table  is  a  brief  summary  of  the  measurements  and  participants. 

All  optical  equipment  was  operated  during  an  intensive  period  from  30  April  to 
6  May  when  the  skies  were  clear  on  every  night  but  one.  Four  aircraft  flights 
were  made  (another  four  flights  had  been  made  earlier  to  calibrate  the  radar) 
and  approximately  35  thermosondes  were  launched  during  this  period.  High  time 
resolution  data  were  logged  on  the  VHF  radar  and  the  sodar  was  operated 
continuously  during  the  optical  measurements.  Synoptic  meteorological 
information  was  carefully  analyzed  and  archived.  A  meeting  was  held  at  AFGL 
in  early  September  to  discuss  the  preliminary  results.  One  of  the  Pi's  (CWF) 

and  three  graduate  students  attended  and  made  presentations. 

2 

The  standard  radar  equation  is  used  to  calculate  Cn  from  the 
received  power.  The  calibration  factor  for  the  radar  is  a  combination  of  a 
number  of  system  constants:  antenna  gain,  line  loss,  receiver  noise,  etc. 
However,  the  antenna  gain  for  the  colinear-coaxial  phased  arrays  used  with  the 
Doppler  radars  has  never  been  rigorously  measured.  Thus,  even  the  most 
careful  determination  of  the  other  calibration  factors  cannot  eliminate  all  of 
the  uncertainty.  Because  of  this,  since  the  Penn  State  radars  are  modeled 
after  the  NOAA/WPL  systems,  we  initially  decided  to  use  the  system  constants 
from  the  NOAA  radars.  The  analysis  of  the  EWAK  data  is  still  ongoing,  but  it 

O 

has  revealed  that  the  radar  calibration  is  quite  good;  the  values  of 
need  to  be  increased  by  about  a  factor  of  two  to  agree  with  the  thermosonde 
balloon  measurements.  A  sample  comparison  of  the  uncorrected  radar, 
thermosonde  and  aircraft  profiles  of  optical  C  is  shown  in  Fig.  2. 


Summary  of  measurements  for  the  EWAK  experiment. 


Measurement 

Institution 

Contact 

Surface  micromet 

PSU 

C. Fairall 

Sodar 

PSU- 

D . Thomson 

VHF1 , VHF2 

PSU 

D. Thomson 

Thermosonde  (CnA2) 

AFGL 

J . Brown 

Aircraft 

ARA 

R.Markson 

Optical  CnA2  profile#l 

RADC 

D .  Stebbins 

Optical  CnA2  profile#2 

AFGL 

E . Murphy 

Optical  scintillometer(rO) 

NPS 

D. Walters 

Optical  isoplanometer 

AFWL 

J .  Davidson 

Analysis  of  the  aircraft  data  has  been  completed  and  published  initially 

in  December  (Beecher,  1987)  in  the  form  of  an  M.S.  thesis.  An  initial 

scientific  publication  has  also  been  completed  (Fairall  et  al. ,  1988).  We 

will  present  a  few  examples  from  this  work.  Fig.  3  shows  a  sample  aircraft 

profile  of  microturbulence  data  from  one  of  the  eight  flights  made  during 

2 

EWAK.  Fig.  4  shows  two  comparisons  of  optical  Cn  obtained  from  the 
aircraft,  the  thermosonde,  the  radar,  the  Van  Zandt  model  and  optical 
scintillometer  data.  The  optical  data  is  highly  variable  over  the  one  hour 
period  but  is  in  agreement  with  the  thermosonde  and  clearly  lower  than  the 
other  two  independent  methods.  Finally,  a  conglomerate  analysis  of  the  small 
scale  turbulence  from  all  flights  is  shown  in  Fig.  5.  Here  the  relation  of 
thermal  and  velocity  microstructure  (as  in  equation  13)  is  examined.  Notice 
that  for  the  special  case  of  temperature 


2  2 

By  plotting  C,j.  versus  (T/g)  36/3 z  Cu  in  log- log  format,  we  expect 
the  points  to  define  the  dependence  of  If  is  constant, 

then  we  expect  a  straight  line  with  a  slope  of  one.  Clearly  this  is  not  what 
we  observed.  For  strong  turbulence  (clearly  active  and  nondecaying  regions 
with  £>10"^) ,  the  data  gives  a  value  of  of  about  0.5,  which  is 
consistent  with  the  jet  stream  value  0.4  of  Kennedy  and  Shapiro  (1980)  and  the 
oceanic  value  0.3  of  Gregg  (1987).  For  the  weaker  turbulence  cases,  this 
value  increases  by  about  a  factor  of  10.  The  significance  of  this  is  not 
clear.  Presently  we  are  not  sure  of  the  effects  of  the  vertical  averaging 
necessary  to  produce  the  300  m  resolution  which  is  required  in  order  for  the 
data  to  be  compatible  with  that  from  the  radar. 

The  second  part  of  the  EWAK  study  is  an  in  depth  investigation  of  the 
2 

radar  and  wind  data  for  the  EWAK  period  by  Mike  Moss,  an  Air  Force 
graduate  student  enrolled  at  Florida  State  University.  This  work  involves  the 
Van  Zandt  model,  the  thermosondes,  and  the  optical  data.  Since  this  is  a 
Ph.D.  thesis  project  (which  requires  three  or  four  years  of  work),  we  do  not 
anticipate  final  results  until  next  year.  A  preliminary  description  of  the 
work  in  progress  has  been  prepared  (Moss,  1986). 

C.  Climatological/Synoptic  Studies 

The  initial  steps  in  the  climatological  study  of  the  synoptic  context  of 

2 

have  been  completed.  We  have  been  routinely  archiving  winds  and  power 
at  one  hour  intervals.  During  a  one  year  period,  each  hour  was  assigned  a 
synoptic  classification  based  on  the  weather  charts.  A  grand  ensemble 
analysis  of  this  data  base  has  not  been  started.  For  the  purposes  of  this 
project,  we  have  completed  a  detailed  analysis  of  one  selected  synoptic 
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feature  that  is  of  particular  interest:  the  jet  stream  and  the  wind  shear  in 

its  vicinity.  This  work  has  now  been  published  initially  in  the  form  of  an 

M.S.  thesis  (Syrett,  1987).  A  scientific  publication  is  now  in  preparation 

(Syrett  and  Thomson,  1988) .  Another  phase  of  this  work  involves  comparison  of 

the  radar  measured  winds  with  both  analyzed  and  predicted  winds  for  the  NWS 

Nested  Grid  Model  (NGM) .  A  paper  on  this  topic  has  been  presented  (Thomson  et 

2 

al.,  1988).  This  work  is  relevant  to  the  numerical  prediction  of  Cn  from 
NGM  products. 

Hourly  measurements  of  wind  speed  and  direction  were  examined  for  two 
prolonged  jet  stream  occurrences  over  western  Pennsylvania.  Data  from  two  of 
the  Penn  State  radars  (McElevys  Fort  and  Crown)  were  stratified  into 
categories  based  on  location  of  the  jet  axis  relative  to  the  site.  Low 
resolution  data  from  the  Crown  radar  were  also  compared  to  the  Pittsburgh 
rawinsonde.  Potential  temperature  profiles  were  obtained  using  isentropically 
interpolated  T  and  Td  soundings .  The  combination  of  measured  wind  and 
interpolated  temperature  profiles  allowed  low  resolution  Ri  profiles  to  be 
generated  for  the  profiler  sounding  volume.  Both  Ri  and  wind  shear  statistics 
were  examined  along  with  pilot  reports  of  turbulence  in  the  vicinity  of  the 
profiler. 

Two  cases  were  examined.  Case  1  lasted  from  7  November  to  14  November 
1986.  Jet  stream  case  2  lasted  from  January  15  to  January  23,  1987.  A  sample 
time-height  cross  section  of  wind  speed  (Fig.  6)  illustrates  the  structure  and 
variability  of  the  jet  stream  at  the  Crown  site  (just  north  of  Cook  Forest  in 
western  Pennsylvania).  The  data  were  separated  into  five  position  categories. 
The  average  wind  speed  and  Ri  profiles  for  case  2  are  shown  in  Fig.  7.  The 
vertical  resolution  of  the  profiler  data  is  slightly  less  than  1  km,  but  was 
interpolated  to  250  m  resolution  for  the  analysis.  The  potential  temperature 
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profiles  also  have  low  horizontal  resolution  because  they  are  extracted  from 
smoothed  fields  from  the  NWS  observing  network.  Notice  the  minima  in  Ri  in 
the  regions  of  maximum  shear  above  and  below  the  peak  of  the  jet.  Pilot 
reports  of  the  altitude  and  severity  of  turbulence  for  the  western 
Pennsylvania  region  were  compiled  for  the  period  of  the  event  (Fig.  8).  The 
PIREPS  were  most  closely  correlated  with  the  magnitude  of  the  shear  observed 
by  the  radar,  rather  than  the  occurrences  of  low  values  of  Ri  (Fig.  9).  This 
figure  unequivocally  shows  that  regions  of  strong  shear  are  almost  sure  to 
have  severe  turbulence.  This  is  consistent  with  (19a)  which  implies  the 
strongest  turbulence  in  the  thickest  layers  (largest  Lq)  and  the  regions  of 
greatest  wind  shear  (8U/<3z) .  Also,  recall  that  (29)  implies  that  the  thickest 
layers  are  also  the  longest  lived  and,  therefore,  the  most  likely  to  be 
encountered  by  aircraft.  This  result  suggests  that,  as  far  as  severe 
turbulence  is  concerned,  it  is  probably  a  mistake  to  attempt  to  derive  the 
Richardson  number  from  low  resolution  sources,  partly  because  regions  of 
strong  shear  are  likely  to  have  small  Ri  anyway. 

Preliminary  results  from  the  comparisons  of  the  profiler  winds  with  the 
NGM  winds  have  been  surprising.  The  profiler  winds  have  been  compared  with 
analyzed  (0  hour  'forecast')  and  predicted  (48  hour  forecast)  NGM  winds.  The 
NGM  data  is  interpolated  to  the  profiler  location  (both  Crown  and  McElevys 
Fort  sites  were  used) .  The  profiler  data  has  been  examined  with  several  time 
resolutions:  the  basic  one  hour  average  at  the  time  of  the  NGM  data,  a  seven 
hour  average  (plus  and  minus  three  hours  about  the  time  of  the  observation) , 
and  a  thirteen  hour  average.  The  comparison  is  best  for  the  longer  time 
average,  indicating  the  amount  of  spatial  averaging  inherent  in  the  regional 
scale  model  or  smoothing  of  extreme  observations  in  the  radar  data.  A  sample 
comparison  for  the  500  mb  wind  speed  and  direction  for  May,  1987,  is  shown  for 


the  analysis  (Fig.  10)  and  the  48  hour  forecast  (Fig.  11).  The  initial 
analysis  is  typically  2.5  m/s  different  from  the  radar  data  and  there  appears 
to  be  no  significant  systematic  bias  in  the  wind  direction.  It  is  believed 
that  the  NGM  underestimates  the  winds  compared  to  the  radar  because  of 
smoothing  of  the  fields  in  the  analysis  procedure.  The  48  hour  forecast  is 
noticeably  more  inaccurate  than  the  analysis,  but  the  deterioration  is 
surprisingly  small. 
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Suggested  Future  Research 


The  general  goal  of  this  research  is  to  investigate  the  fundamental 

relationsnips  between  mean  atmospheric  structure  and  microturbulence 

2 

parameters  such  as  Cn  .  We  wish  to  study  the  processes  that  lead  to 

2 

temporal/horizontal  variability  and  vertical  distribution  of  Cn  .  The 
next  obvious  step  in  this  investigation  is  to  acquire  and  analyze  data  from 
three  primary  sources : 

(1)  Existing  data.  Radar  data  from  VHF1  and  VHF2  archived  over  the  past 
two  years  will  be  available  for  completion  of  a  climatological 
analysis.  Raw  aircraft  turbulence  data,  which  is  archived  on  FM  tapes, 
and  raw  aircraft  mean  profile  data,  which  is  archived  on  digital  tapes 
with  2  second  resolution,  from  the  EWAK  experiment  could  be  reanalyzed 
with  finer  resolution. 

(2)  New  Profiler  data.  We  will  continue  to  archive  one  hour  average 
radar  data  from  VHF1 ,  VHF2 ,  and  VHF3  (which  is  to  be  in  operation  in 
September,  1988).  We  will  also  archive  one  hour  average  data  from  the 
UHF  radar,  which  is  to  begin  regular  operation  at  State  College  in 
July,  1988.  Also,  we  hope  to  begin  operation  of  the  mm-wave  radiometer 
system  in  State  College  about  July  of  1988.  This  system  will  archive 
two  minute  time  resolution  measurements  of  integrated  water 
vapor/liquid  and  temperature  profiles.  The  mesoscale  triangle  (about 
140  km  sides)  of  VHF  radar  combined  with  the  UHF  and  thermodynamic 
sounders  will  provide  an  unprecedented  look  at  atmospheric  turbulence 
and  dynamical  processes. 

(3)  Intensive  Turbulence  and  Optical  Propagation  Experiment  (EWAK  II). 


Sometime  in  the  near  future  it  makes  sense  to  have  another  field 
program  similar  to  EWAK.  We  would  try  to  work  with  the  same 
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institutions  as  before  for  the  optical  measurements.  We  could  also  try 
to  interest  and  invite  NASA  and  NOAA  groups  with  relevant  measurement 
systems  (e.g.,  Melfi's  raman  lidar) . 

The  variety  of  scientific  issues  of  interest  were  discussed  in  section 
III.  We  can  divide  the  issues  into  four  general  categories:  spatial 
structure,  turbulence  dynamics,  climatology/synoptic  correlations,  and 
comparisons  of  optical,  radar  and  in  situ  measurements. 

A.  Spatial  Structure 

The  vertical  distribution  of  the  active  turbulence  areas  can  be  studied 

with  combined  aircraft  and  radar  data.  We  are  interested  in  the  probability 

distributions  of  patch  size,  patch  lifetime,  turbulence  scales, and  the 

autocorrelation  functions  of  patch  distribution.  We  can  also  use  aircraft 

measurements  to  examine  the  probability  distributions  of  lapse  rate  and  shear 

(basic  parameterizations  in  the  Van  Zandt  model) .  The  horizontal  and  temporal 

persistence  of  individual  turbulent  layers  can  be  studied  by  cross  correlation 
2 

the  Cn  data  from  the  three  VHF  radars  and  by  examining  level  flight 
aircraft  data. 

B.  Turbulence  Dynamics 

Here  the  emphasis  is  on  the  microturbulence  scaling  relations  (e.g.,  the 

2  2 

ratio  C.j,  /C  ),  the  generation  and  decay  process,  and  the  turbulence 

size  scales  (L^  and  Lq) .  Aircraft,  radar  and  thermosonde  data  can  all  be 

used.  We  suggest  investigating  the  use  of  the  radar  to  measure  e  from  the 

width  of  the  Doppler  spectrum.  Furthermore,  we  can  compare  e  derived  from  the 

width  of,  for  example,  a  300  m  range  resolution  spectrum  with  that  from  the 

mean  shear  derived  from  overlapping  100  m  range  resolution  data.  Thus,  we  can 

2  2 

examine  simultaneous  measurements  of  and  .  The  size  scales  of 
turbulence  can  be  investigated  in  several  ways.  The  operation  of  colocated 
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VHF  and  UHF  radar  will  also  provide  information  about  the  relative  intensity 
of  turbulence  size  scales  since  they  scatter  from  different  atmospheric  scales 
(A/2=«3  m  and  0.3  m  respectively).  Profiles  of  mean  temperature  from 
aircraft,  rawinsonde,  thermosonde,  and  radiometers  would  give  an  indication  of 
N  which,  combined  with  £, gives  L^.  The  vertical  velocity  spectrum  from  the 
radar  can  also  give  an  indication  of  size  scales.  In  active  turbulent  layers 
it  should  reveal  Lo;  in  more  quiescent  layers  it  should  indicate  N.  In  the 
intensive  experiment,  a  great  deal  of  effort  could  be  devoted  to  examining 
these  issues. 

C.  Climatology/Synoptic 

The  mesoscale  triangle  of  VHF  radar  will  provide  high  quality  dynamical 

variables  usually  not  available  (e.g,  divergence  and  vorticity)  that  will  be 

2 

extremely  useful  in  interpreting  variability  in  Cn  .  For  example,  the  jet 
stream  study  referred  to  in  section  IV  would  have  greatly  benefited  from  a 
more  quantitative  classification  of  conditions  than  the  crude  'distance  from 
the  jet  axis'  categories  used.  The  temperature/humidity  radiometers  also 
provide  very  accurate  est^Jls^gigofdYnami  -  height  (e.g.,  500  mb)  to  augment 
the  radar  information.  We  presently  plan  to  continue  the  large  scale  synoptic 
classifications  for  the  gross  climatological  analysis.  We  will  also  continue 
the  studies  of  the  comparisons  of  radar  data  and  NGM  analyzed  and  forecast 
fields . 

D.  Comparison  of  Optical  and  Meteorological  Sensors 

We  strongly  recommend  that  another  joint  meteorological  and  optical 
experiment  be  held  a  Penn  State.  EWAK  I  revealed  that  there  are  still 
substantial  differences  between  the  various  observing  systems  (often  an  order 
of  magnitude) .  Some  of  this  is  due  to  differences  in  temporal  and  spatial 
averaging  techniques  but  a  great  deal  of  the  disagreement  must  be  due  to 
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undiagnosed  problems  with  the  sensors  and  data  processing.  A  new  intensive 
program  would  also  provide  us  with  other  sources  of  data  (eg.,  thermosondes 
and  aircraft)  and  much  better  spatial  and  temporal  resolution  of  the  important 
variables.  This  is  necessary,  in  particular,  to  determine  the  importance  of 
the  averaging  process  on  the  interpretation  of  measurements  and  models. 
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Figure  5a.  A  comparison  of  simultaneous  values  of  temperature  and  velocity 
microturbulence  when  normalized  using  equation  13.  A  slope  of  1  would 
indicate  that  was  a  constant. 
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Figure  5b.  A  composite  comparison  of  normalized  by  Cu  and 
potential  temperature  gradient  (as  in  Eqs .  13  and  30)  as  a  function  of  the 
activity  parameter,  A.  A  constant  value  would  indicate  that  7 g  was  a 
constant . 
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Figure  la.  Average  wind  speed  profiles  with  the  measurement  site  within  100 
km  of  the  jet  axis.  The  horizontal  bars  indicate  the  standard  deviation  for 
the  period.  The  upper  panel  is  case  1;  the  lower  panel  is  case  2. 
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Figure  9.  Surface  pLots  of  wind  shear  above  the  Crown  radar  during  21  January 
1987.  The  surface  defines  the  time,  altitude  and  magnitude  of  the  wind  shear. 
The  shaded  areas  indicate  the  pilot  reports  of  turbulence:  dark  implies  severe 
and  dotted  implies  moderate.  The  bottom  panel  is  a  90-degree  rotation  of  the 
upper  panel  to  provide  a  different  view.  Note  that  all  pilot  reports  or 
severe  turbulence  occurred  during  periods  of  maximum  wind  shear. 
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APPENDIX  A 


Profiler  System  Description 

Mesoscale  meteorological  measurements,  analysis  and  prediction  are  some 
of  the  principal  areas  of  research  in  the  Department  of  Meteorology  at  Penn 
State.  For  more  than  a  decade  those  members  of  the  faculty  concerned  with 
mesoscale  analysis,  numerical  modeling  and  forecasting  have  been  frustrated  by 
the  spatial  and  temporal  inadequacy  of  conventional  network  observations  for 
both  research  and  operational  applications.  For  more  than  five  years  the 
Department  had  sought  the  substantial  financial  resources  required  to  deploy  a 
network  of  VHF  Doppler  (ST)  radars  and  millimeter  wave  radiometers  for 
"operational"  test  and  evaluation  for  wind  and  thermodynamic  profiling. 
Construction  of  the  ST  radar  network  began  in  fall  of  1983  using  funding 
provided  by  the  Air  Force  Office  of  Scientific  Research  (through  the  DoD 
University  Research  Instrumentation  Program)  and  the  University.  In  1986  a 
second  AFOSR  URIP  grant  was  obtained  to  fund  the  acquisition  of  temperature 
(50-60  GHz)  and  humidity  (20-30  GHz)  profiling  radiometers.  This  document 
will  only  discuss  the  radar  systems. 

For  the  foreseeable  future  the  Penn  State  ST  radar  program  will  be 
focused  on  applications  rather  than  systems  development  research.  Deployment 
of  the  systems  would  not  have  been  possible  without  the  outstanding 
cooperation  provided  by  C.G.  Little  and  R.  Strauch  and  their  colleagues  at  the 
Wave  Propagation  Laboratory,  and  also  J.  Brosnahan  of  Tycho  Technology  from 
whom  we  bought  all  of  the  receivers,  transmitters  and  antennas.  With  regard 
to  the  other  major  systems  components,  we  have  assembled  in-house,  from  WPL 
documentation,  the  time -domain- integrator  and  computer  interfaces  and  have 
purchased  WPL  software-compatible  Data  General  Corp.  Eclipse  computers  for 
each  system. 
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The  Penn  State  network  consists  of  three  6  m  wavelength (VHF)  and  one  0.7 
m  wave 1 eng th(UHF)  radars.  Fig.  1  indicates  the  approximate  location  of  a 
mesoscale  triangle  formed  by  the  three  VHF  radar  sites  within  the  routine 
rawinsonde  network.  In  Table  1,  a  short  summary  of  the  specifications  for  the 
four  radars  is  given.  Experiences,  plans  and  improvements  for  the  PSU  network 
are  summarized  below. 

A.  VHF1  50-MHz  radar  located  15  km  south  of  State  College,  PA. 

1.  This  system  became  fully  operation  June  27,  1985.  The  primary  reason 
for  system  failures  since  the  onset  has  been  AC  outages  which  are  prevelant  in 
this  area.  Battery  back-up  and  computer-controlled  autorestart  of  the 
transmitters  has  circumvented  this  problem. 

2.  Initial  performance  statistics  done  by  Frisch  et  al.,  WPL  on  August 
1985  data  ,  indicate  very  good  performance  by  VHF1.  On  Beam  #1  the  next  to 
last  range  gate  (16.8  km  MSL)  was  able  to  make  a  wind  measurement  99%  of  the 
time,  while  Beam  #2  was  able  to  measure  the  wind  at  this  height  85%  of  the 
time.  The  difference  befrce"  beams  is  probably  due  to  better  pickup  by 
beam  #2  of  computer/electronic  noise  from  the  building.  Samples  of  data  from 
the  low  altitude  range  of  VHF1  are  shown  in  Figs.  2  and  3. 

3.  A  vertical  beam  was  added  in  March,  1987,  and  the  latest  NOAA/WPL 
software  was  installed. 

B.  VHF2  50-MHz  radar  located  in  NW  Pennsylvania  near  Crown,  PA. 

1.  This  system  was  installed  May  1,  1986  and  has  operated  since. 
Performance  is  often  slightly  better  than  VHF1  but  a  source  of  interference 
(radiotelephones  from  a  local  trucking  company)  leads  to  periods  of  reduced 


data  quality. 


2.  A  vertical  beam  was  added  in  January,  1987,  and  the  latest  NOAA/WPL 
software  was  installed. 

C.  VHF3  50-MHz  radar  located  in  SW  Pennsylvania  near  Somerset,  PA. 

1.  All  hardware  is  ready  for  installation  pending  selection  and 
preparation  of  a  site. 

2.  A  three  beam  system  similar  to  VHF1  and  VHF2  will  be  used.  Some  of 
the  hardware  from  Tycho  Technology  has  been  slightly  updated. 

D.  UHF1  portable  405 -MHz  radar  to  be  semi -permanently  based  at  the  PSU 
Circleville  Farm  in  downtown  State  College,  PA. 

1.  This  system  has  been  used  on  three  major  field  deployments:  the 
SPACE/MIST-COHMEX  experiment  in  Alabama  (June-July,  1986),  the  FIRE 
stratocumulus  IFO  on  San  Nicolas  Island,  CA  (July,  1987),  and  the  Arizona 
monsoon  experiment  (August,  1987). 

2.  Preparations  for  basing  at  Circleville  Farm  are  to  be  completed  in 


July,  1988. 


Table  1:  Specifications  for  the  Penn  State  ST  Radars 


Item 

VHF1.VHF2.VHF3 

UHF1 

Type 

Pulsed  Doppler,  3  beam 

Pulsed  Doppler,  3  beam 

Location 

1:  S.  State  College 

2:  Crown,  PA 

3:  Somerset, PA 

In  State  College,  trans 
portable 

Frequency  (nom.) 

49.9  MHz 

404  MHz 

Bandwidth 

300,  100  KHz 

1,  0.3,0.07  MHz 

Peak  Power 

30  KW 

30  KW 

Pulsewidth 

3.67,  9.67  fise c 

1,  4,  16  /xsec 

Antenna : 

Type 

Dimensions 

Zenith  angles 

Phased  array  CoCo 

50  m  X  50  m 

90* (vert) ,75° 

Phased  array  CoCo 

9  m  X  9  m 

90° (vert),  75“ 

One  site  computer 

Data  Gen.  Eclipse 

Data  Gen.  Eclipse 

On  site  processing* 
at  PW-  : 

3.67  /is  9.67  /is 

l//s  4/xs  16/is 

Time  domain  ave 
Spectral  ave 

Height  spacing  (m) 
Spec,  resolut. (m/s) 
Maximum  absolute  rad. 
velocity  (m/s) 

400  125 

8  16 

290  870 

0.49  0,31 

15.7  19.6 

112  70  35 

16  32  64 

100  300  800 

0.29  0.29  0.29 

18.25  18.25  18.25 

*Software  controllable 
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NOTES  CONCERNING  THE  USE  OF  CLEAR  AIR  DOPPLER  RADARS 


FOR  THE  MEASUREMENT  OF  Cu 2 


Robert  M.  Peters 


Department  of  Meteorology 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Recently  the  Penn  State  clear  air  radar  network  has  beer, 
increasingly  utilized  not  only  for  wind  profiling,  but  also  for  the 
measurement  of  the  refractive  structure  coefficient  associated  with 
isotropic  turbulence,  Cu 2  .  This  coefficient  is  related  to  received 
power  amplitude  of  active  indirect  sensing  systems  when  signal 
backscatter  is  caused  solely  by  refractive  index  variations  due  to 
isotropic  turbulence.  The  physical  mechanisms  associated  with  this 
backscatter  in  the  atmosphere  are  described  in  detail  by  WESELY 


& 


Vi 


(1976)  . 

The  present  Penn  State  clear  air  radar  network  consists  of  two 
VHF  ( 5M  wavelength)  Doppler  radars  and  a  transportable  UHF  (75  CM 
wavelength)  Doppler  radar.  The  systems  are  designed  primarily  for 
wind  profiling  of  the  troposphere  and  stratosphere.  V"  a  Penn  State 
network  is  described  by  THOMSON  et  al.  (1934)  and  is  based  upon  ohe 
concepts,  signal  processing  hardware  designs  and  software  provided  by 
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the  NOAA/ERL  Wave  Propagation  Lab  (STRAUCH  et  al . ,  1934).  An  ongoing 
experimental  research  effort  at  Penn  State  pertains  to  measurements  of 
turbulent  backscatter  from  radars,  optical  systems,  acoustic  systems 
and  direct  measurements  of  turbulence  via  aircraft  and  thermosondes. 
During  a  recent  experiment,  one  VHF  radar  was  used  for  this  purpose 
and  measurements  of  the  VHF  structure  coefficient  were  obtained  (MOSS, 
1986) .  It  has  become  evident  that  the  quantitative  measurement  of 
radar  system  power  parameters  that  are  required  for  the  determination 
of  Cn 2  are  prone  to  errors  that  are  frequently  not  addressed  in 
designs  that  are  optimized  for  Doppler  measurements. 


In  order  for  backscattered  power  to  be  related  to  isotropic 
turbulence,  it  must  be  known  that  the  received  signal  amplitude  is  not 
modified  by  other  backscatter  sources,  reflections  or  interference 
from  other  emitters.  Continuous  interference  is  usually  man  made  from 
sources  such  as  arcing  power  lines,  radio  transmissions  or  broadband 
noise  from  the  digital  systems  of  the  radar.  Good  site  selection  and 
radar  system  engineering  practices  will  prevent  interference  from 
these  types  of  emission.  Hard  targets  will  cause  undesired 
reflections.  Reflections  due  to  ground  clutter  may  be  partially 
eliminated  in  spectral  processing  but  are  best  avoided  by  choosing 
radar  sites  that  have  no  large  targets  visible  at  distances  comparable 
to  the  radar's  range  response.  Hard  reflections  from  aircraft  will 


override  the  atmospheric  backscattered  signal 


The  consensus 


algorithm  used  for  averaging  in  the  Penn  State  radar  software  (STRAUCH 
et  al . ,  1984)  will  eliminate  the  short  term  occurrence  of  aircraft 
reflections  from  an  averaged  data  .'.et.  However,  the  location  cf  major 
airways  should  be  a  factor  during  radar  site  selection.  Enhanced 


signal  returns  from  the  atmosphere  from  other  natural  sources  of 


anisotropic  irregularities  will  occur  with  a  zenith  pointing  beam  as 
described  by  DOVIAK  and  ZRNIC '  {1983).  As  shown  by  DOVIAK  and  ZRNIC’ 
(1983) ,  the  response  of  VHF  radars  to  anisotropic  irregularities  will 
be  insignificant  at  a  15  degree  off  zenith  pointing  angle,  assuming 
that  the  layered  structure  causing  the  signal  enhancement  is 
horizontal.  At  Penn  State,  one  of  the  15  degree  off-axis  beams  of  the 
VHF  radar  antenna  is  used  for  the  measurement  of  Cn 2  as  well  as  for 
one  wind  velocity  component.  Antenna  beam  side  lobe  response  to 
enhanced  signal  returns  from  the  vertical  are  further  reduced  by 
having  a  zenith  null  designed  into  the  antenna  pattern. 

In  order  to  make  quantitative  measurements  of  received  signal 
amplitudes,  the  radar  systems's  hardware  must  be  calibrated  or  its 
calibration  inferred  from  other  sources.  System  hardware  for 
calibration  purposes  are  the  antenna  system,  transmitters,  receivers 
and  analog  to  digital  converters  used  for  parameter  measurements. 

The  antenna  patterns  of  large  upward  pointing  VHF  phased 
arrays,  such  as  those  used  at  Penn  State,  are  difficult  to  measure  in 
their  far  fields.  Computer  antenna  pattern  models  are  quite  accurate 
with  regard  to  antenna  pattern  shape.  Experience  has  shown  that  an 
antenna  will  usually  have  a  gain  within  a  few  decibels  (dB)  of  that 
predicted  by  a  validated  antenna  model,  if  the  antenna  is  well 
engineered  and  installed.  Pattern  altering  metallic  objects  should 
not  be  placed  in  or  near  the  antenna  field  if  reliance  is  placed  upon 
previously  determined  antenna  gain  figures  that  are  to  be  used  in  the 
measurement  of  structure  parameters.  Astronomical  noise  sources  have 
been  used  to  confirm  antenna  pointing  maxima.  Galactic  sources  are 
useful  in  the  northern  hemisphere  (MOSS,  1986)  and  the  sun  can  be  used 


in  equatorial  regions  (RIDDLE,  1985).  Precise  patterns  and  gains  have 
been  determined  for  the  Mu  radar  facility  in  Japan  by  integrating 
direct  satellite  measurments  as  described  by  SATO  et  al  (1985).  The 
phase  and  power  distribution  provided  to  the  elements  of  the  VHP 
phased  array  may  be  directly  measured  to  provide  a  more  accurate  input 
to  modeled  antenna  patterns.  The  performance  of  an  uncalibrated 
antenna  system  may  be  statistically  compared  with  other  systems  of 
known  performance  and  its  gain  inferred  (MOSS,  1986),  as  long  as  other 
system  components  are  calibrated.  Computer  estimation  of  Ct!  2  may  be 
derived  from  atmospheric  soundings  and  used  to  estimate  system, 
calibration.  A  field  tested  computer  model  developed  by  VAMZANDT  et 
al .  (1981)  could  possibly  be  used  for  this  purpose. 

Transmitter  parameters  such  as  power  and  pulse  width  must  be 
known.  The  Penn  State  network  incorporates  transmitters  which  have 
internal  microprocessor  control  and  monitoring  systems.  The  radar 
system  executive  control  computer  at  each  site  provides  a  means  to 
remotely  monitor  and  control  each  radar  transmitter  in  the  network. 
Transmitter  power  may  be  monitored  via  this  arrangement,  however,  this 
indirect  method  of  measurement  is  prone  to  error  if  the  associated 
electronic  circuits  are  not  calibrated.  With  Penn  State  systems,  the 
transmitted  power  level  is  derived  from  an  adjustable  RF  voltage 
monitoring  circuit  with  its  output  supplied  to  the  analog  to  digital 
converter  of  each  transmitter's  processor.  Each  transmitter  includes 
an  integral  directional  coupler.  To  ensure  the  accuracy  of  the 
transmitter  power  monitor,  the  coupling  coefficient  of  each 
directional  coupler  should  be  calibrated  and  power  measurements  taken 
to  check  and  calibrate  the  microprocessor  based  power  monitoring 


The  receiver  system  includes  not  only  the  radar  receivers  and 
preamplifiers,  but  also  the  analog  to  digital  converters.  For 
structure  parameter  measurements,  the  receiving  system  gain  must  be 
calibrated.  For  relative  power  measurements  of  received  signal 
amplitude,  the  calibration  must  be  done  with  test  instrumentation. 
Any  changes  to  the  receiving  system,  such  as  replacement  of  a  receiver 
preamplifier,  would  require  recalibration.  For  radar  systems  such  as 
the  Penn  State  UHF  Doppler  radar,  this  is  the  most  practical  method  of 
calibration  since  the  UHF  system  utilizes  one  receiver  time  shared 
among  several  antenna  pointing  angles.  The  minimun  discernable 
received  signal  power  of  the  UHF  system  is  limited  by  the  100  K  noise 
temperature  (approximate)  of  the  receiver  preamplifier.  This  type  of 
calibration  is  unnecessary  for  the  VHF  Doppler  radars.  For  these 
systems,  the  minimum  signal  detectabilty  is  limited  by  a  galactic 
noise  temperature  backround  in  the  order  of  several  thousand  Kelvins 
(KO,  1958).  For  this  reason,  the  galactic  noise  level  is  used  as  an 
automatic  calibration  for  all  receivers  at  Penn  State  VHF  radar  sites. 
This  is  accomplished  by  archiving  a  spectrally  derived  estimate  of 
consensus  averaged  signal  to  noise  ratio  (most  systems  presently 
archive  relative  power  on  a  routine  basis) .  Structure  parameter 
errors  at  close  ranges  using  VHF  radar  have  been  observe  during  recent 
experiments  (MOSS,  1986).  The  cause  for  this  class  of  error  has  not 
been  studied  but  appears  either  to  be  due  to  a  near  field  antenna 
response  error  or  receiver  dynamic  range  limiting. 

Less  than  ideal  receiver  low  pass  filter  stability, 
calibration  and  adjustment  has  been  the  most  significant  cause  of 
error  of  spectrally  derived  structure  parameters  from  the  VHF  radars  . 
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Specifically,  the  presence  cf  DC  offsets  and  gain  imbalance  be*v-?n 
real  and  quadrature  channels  will  produce  a  host  of  errors  in  spectral 


signal  processing, 


Each  received  signal  spectrum  is  actually  the 


result  of  several  averaged  spectra,  each  of  which  is  derived  from  a 
time  series  of  coherently  integrated  received  signal  amplitudes  at 
each  radar  altitude.  Each  signal  spectrum  is  derived  from  the  complex 
Fast  Fourier  Transform  of  the  receiver’s  range  gated  real  and 
quadrature  outputs.  Peaks  that  may  occur  in  each  spectrum  about  the 
zero  frequency  line  at  an  equal  amplitude  at  all  range  gates  are 
usually  the  result  of  the  presence  of  a  mean  or  trend  in  the  complex 
signal  data  sef  that  is  not  completely  removable  by  the  signal 
processing  software.  Image  peaks  that  may  occur  opposite  the  spectral 
peak  associated  with  the  clear  air  signal  return  usually  are  a  result 
of  dynamic  range  limiting  of  the  receiving  system.  A  consistent-  i m  => g 
peak  in  all  range  gates  is  usually  associated  with  a  gai’-  cr 
quadrature  imbalance  between  the  I  and  Q  channels  of  the  receiving 
system.  Note  that  the  receiving  system  also  includes  the  analog  to 


digital  converters  of  the  signal  processing  system. 


Imbalances 


between  the  I  and  Q  analog  to  digital  converters  will  produce  the  sa: 


errors 


Errors  produced  by  means,  trends  and  imbalances  produce 


errors  in  structure  parameter  data  that  are  not  as  easily  removed  as 
compared  to  Doppler  velocity  measurements.  The  amplitude  of  a 
spectral  image  will  always  be  lower  than  the  true  signal  so  that 
Doppler  peak  detecting  algorithms  will  never  choose  the  image  peak  for 
Doppler  processing.  Any  means  or  trends  will  be  have  an  effect  in 
Doppler  velocity  processing  only  if  the  zero  amplitude  peak  is  greater 
than  the  clear  air  peak.  If  the  zero  peak  is  larger,  then  the 


artificially  constant  velocity  produced  in  the  beam  component  (which 
is  frequently  blamed  on  ground  clutter)  can  be  easily  removed  from 
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c  ' 


wind  velocity  data  sets  by  computer  post  processing.  The  effect  of 
these  types  of  spectral  errors  in  obtaining  received  signal  power  or 
signal  to  noise  ratio  are  not  as  obvious.  The  spectral  processing 
algorithm  will  include  ground  clutter  peaks  as  part  of  the  returned 
backscattered  signal  at  low  radial  wind  velocities.  Any  image  peaks 
that  occur  will  also  bias  the  estimates  of  received  signal  power  and 
possibly  produce  a  velocity  bias  to  the  estimate.  The  presence  of 
spectral  images  will  contribute  to  the  estimated  noise  level  thereby 
producing  errors  in  the  signal  to  noise  ratio. 


To  reduce  the  errors  associated  with  the  receiver  system  for 
all  types  of  measurements,  the  receivers  must  be  carefully  adjusted  in 
the  lab  for  gain  balance  and  zero  offsets  after  a  period  of 
stabilization .  Software  solutions  can  contribute  to  the  reduction  of 
error.  The  present  software  used  at  Penn  State  (STRAUCH  et  al . ,  1334) 


removes  the  mean  from  the  received  signal  data  set.  Trend  removal  is 
more  difficult  for  the  real  time  nature  of  our  signal  processing 
system  however  improvements  have  been  suggested  by  STRAUCH  (1536). 
The  ultimate  solution  is  to  treat  the  disease,  not  the  symptoms.  The 
filter  and  output  stages  of  VHP  and  UHF  Doppler  receivers  should  be 
designed  for  maximum  amplitude  stability  as  well  as  frequency 
stability.  The  manufacturer  of  the  radio  frequency  subsystems  of  the 
Penn  State  Doppler  radars  are  now  providing  improved  receivers  for  new 
systems  and  retrofits  for  present  systems  to  reduce  the  possibility  of 


these  types  of  problems 


The  Penn  State  UHF  system  will  have  a 


programmable  attenuator  added  to  the  receiver  system.  This  will  allow 
the  receiver's  gain  to  be  fixed  at  an  ideal  level  depending  upon 


operating  requirements 


This  feature,  along  with  improved  on-line 


spectral  display  software,  will  allow  field  optimization  of  radar 
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parameters  to  reduce  the  sources  of  error  described  herein.  An 
improved  method  of  spectral  signal  detection  has  also  been  developed 
(PETERS  and  WILLIAMS,  1988).  The  new  method  provides  a  better  signal 
peak  selection  mode  in  the  presence  of  ground  clutter  or  in  low  signal 
to  noise  conditions.  Careful  attention  to  clear  air  radar  system 
performance  will  result  in  better  estimates  of  all  measurable  clear 
air  parameters. 
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